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Applying scale-free mass estimators to the Local Group in 
Constrained Local Universe Simulations 
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ABSTRACT 

We use the recently proposed scale-free mass estimators to determine the masses of the 
Milky Way (MW) and Andromeda (M31) galaxy in a dark matter only Constrained Local 
UniversE Simulation (CLUES). While these mass estimators work rather well for isolated 
spherical host systems, we examine here their applicability to a simulated binary system with 
a unique satellite population similar to the observed satellites of MW and M3 1 . We confirm 
that the scale-free estimators work also very well in our simulated Local Group galaxies with 
the right number of satellites which follow the observed radial distribution. In the isotropic 
case and under the assumption that the satellites are tracking the total gravitating mass, the 
power-law index of the radial satellite distribution A^(< r) oc r^^'' is directly related to the 
host's mass profile M{< r) (x r^^" as a = 7 — 2. The use of this relation for any given 7 
leads to highly accurate mass estimations which is a crucial point for observer, since they do 
not know a priori the mass profile of the MW and M3 1 haloes. We discuss possible bias in the 
mass estimators and conclude that the scale-free mass estimators can be satisfactorily applied 
to the real MW and M3 1 system. 
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1 INTRODUCTION 

Although measurements of gas rotation curves are often precise 
enough to constrain the inner most mass of galaxies like the Milky 
Way (MW) and Andromeda (M31) (within a few tens of kpc), kine- 
matics of a tracer populations are needed to compute the mass 
within greater radii. These tracers can either be globula r clusters 
or pla netary nebulae (e.g. Schuberth et al. (2010):Woo dlev et al.l 
l 2010l)). halo stars jXue et al.l |2008|) or satellite galaxies (e.g. 



IWatkins et all 1 I2OIOI) ). Since the kinematics of these objects are 
determined by the underlying host potential they allow for an es- 
timate of the enclosed mass within their respective distances from 
the center of the host. 

Kinematic data of galaxies in clusters have already been used 
to c ompute the mass profil es and galaxy orbits in nearby clus- 
ters ( IWoitak&Lokasll2010h : moreover, the mass of four Milky 



Way dwarf spheroidals (dSphs) satellites were constra ined with 
high precision thanks to kinematic data sets ( lLokasl2009l) . Line-of- 
sight kinematic observations enable accurate mass determinations 
at half-light radiu s for spherical galaxies such as the MW dSphs 
jWolf etalJboiOl) : at both larger and smaller radii however, the 
mass estimation remains uncertain because of the unknown veloc- 
ity anisotropy. 

Regarding our own Galaxy, having position and proper mo- 
tion data of the MW's satellite galaxies would allow one to satis- 
factorily apply the great majority of kinematic mass estimators to 
the calculus of the Milky Way's mass, includ ing the recently pro- 
posed "scale-free projected mass estimator" dWatkins et al.ll201ol 
hereafter WIO). 

In the very near future the knowledge of the full six- 
dimensional phase-space information for all objects, in the close 
Universe, brighter than G ~ 20 mag, is goingto be dramatically 
improved thanks to space missions, like GAIj^ whose goal is to 
create the largest and most precise three dimensional chart of the 
Milky Way by providing precise astrometric data like positions. 
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parallaxes, proper motions and radial velocity measurements for 
about one billion stars in our Galaxy and throughout the Local 
Groi ip (LG). 

lAn et alj ( 1201 2|) recently showed that new proper motions data 
with the targeted GAIA accuracies will be able to outperform the 
presently existing line-of-sight based mass estimators. But until the 
proper motions of these satellite galaxies become available, one 
needs to rely on assumptions and simplifications. 

One of the first estimators of the m ass contained wi t hin th e 
LG is based on the "timing " argument of iKahn & WoltieJ l ll959l) . 
More accurate mass estimators for spherical systems are based ei- 
ther on the virial the orem or on the moments of projected mass, 
as first introduced by ISahcall & Tremain el (ll98lh . They assumed 
that only projected distances and line-of-sight velocity informa- 
tion were available, and demonstrated the goodness of the projected 
mass estimator. The main advantages of such a projected mass es- 
timator over the virial theorem, neglecting the uncertainties in the 
eccentricity distribution, are that they are unbiased, their variance 
is known, and they converge to the real mass with an error pro- 
portional to N^^'"^, where TV is the sample number. Moreover, the 
information from every tracer particle is equally weighted, contrary 
to what happens for the virial theorem case. 

Previous studies successfully used these mass estimator meth- 
ods to "weigh" M31; and more recently, WIO developed alterna- 
tive forms of estimators that can also be applied to the calculus 
of the MW's mass: they rely on the assumption that both the host 
galaxy and its distribution of tracer objects are spherically symmet- 
ric. What is still unclear however, is the shape of the MW and M31 
halo, with various authors in the literature disagreeing over whether 
its tr iaxial (Law et al. 200 ^ or spherical IKopqsov et al. 2 010l) . 

iDeason et al] kOllh and lEvans et alj j201lh have demon- 
strated the statistical validity of the WIO mass estimators using 
a set of 431 parent haloes and 48 64 associated sate llite galaxies, 
taken from the GIMIC simulations l lCrain et al .l2009h : under the as- 
sumption of having a host profile of the type NFW dNavarro et al.l 
Il996h . they found that the fraction of estimated halo mass which 
lies within a factor of two of the tme mass is about 80%. 

In this work we aim to gauge the quality of the method 
introduced in WIO by using the Local Group identified in the 
WMAP5 dark matter only constrained cosmological simulation of 
the CLUES project, a numerical laboratory for testing the appli- 
cability of such a method to the MW and M3 1 . 

Observational data of the nearby Universe are used to con- 
strain the initial conditions of the CLUES simulations. These con- 
strained simulations, in which the Local Group lies in the right cos- 
mological environment, provide a complementary approach, with 
respect to cosmological simulations, to make a comparison be- 
tween numerical results and observations. Thus, verifying the ro- 
bustness of the WIO mass estimators in our unique simulated LG is 
an important test in addition t o the more statistica l methods offered 
by cosmological simulations dOeason et al.ll201lh . 

The idea is to verify whether these estimators can accurately 
be applied to a system such as the one found in our LG and 
composed of the Milky Way and the Andromeda galaxy. The 
arrangement and formation history of this galactic binary sys- 
tem, according to our present state-of-the-art of numerical simu- 
lations, is rather uni que and involves pref erential infall directions 
of t heir subhaloe s ( Libeskind et al.l I2OI ll) , a backsp lash popula- 
tion dKnebe et ai]|201lh . and even renegade satellites dKnebe et al.l 
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I20T ij), i.e. satellites that change their affiliation from one of the two 
hosts to the other. Furthermore, the MW and M31 satellites do in 
fact re member the non-rand om nature of their infall after several 
orbits dLibeskind et al.l2012l) . We also need to mention that - when 
comparing constrained against un-constrained simulations - only 
1-3% of the Local Group candida tes share similar formation prop- 
erties ( iForero-Romero et al .11201 ih . 

Moreover, the observed Milky Way satellites are found to be 
highly anisotropical, lying within a thin disc which is inclined 
with respect to the MW's one, with a minor-to-major axis ratio 
c/a ~ 0.3: this flattened distribution is not compatible with the 
sate llites to have been randomly selected from an isotropic sub- 
set dKroupa et allbOOSl i lMetz et alj|2007l , l2008l) . Previous cosmo- 
logical simulations showed anisotropy in the subhaloes popula- 
tion, with the brightest satellites distributed along dis k-like struc- 
tures, consistently with the obser ved MW satellite s (cf. lKnebe et al.l 
12004: Libeskind et al, 2005; Z entner et al.ll2005l) . This anisotropy, 
which is also observed in our simulated subhaloes, may in princi- 
ple cause a bias in the application of the mass estimator, since the 
hypothesis of spherical symmetry is broken. 

We therefore raise (and answer) the question about the appli- 
cability of scale-free mass estimators to such a special system as 
the Local Group. 



2 THE SIMULATION 

Here we briefly describe the simulation and the way (sub-)haloes 
have been identified within it. 



2.1 Constrained Simulations of the Local Group 

The dark matter only simulation used here forms part of the Con- 
strained Local UniversE Simulations (CLUES) project and is de- 
signed to re produce the Local Gr oup of galaxies within a WMAP5 
cosmology dKomatsu et ai]| 20091') , i.e. = 0.279, Q.^ = 0.046, 
Oa = 0.721. We use a normalization of cts = 0.817 and a slope 
of the power spe ctrum of n = .96. We used the TreePM A^'-body 
code GAD GET 2 dSpringefcoOSi) to simulate the evolution of a cos- 
mological box with side length of Lbox ~ 64/i^^Mpc. Within 
this box we identified (in a lower-resolution run utilizing 1024'^ 
particles) the position of a mo del local group that cl osely resem- 
bles the real Local Group (cf. iLibeskind et al.ll2010h . This Local 
Group has then been re-sampled with 64 times higher mass reso- 
lution in a region of 2/i^^Mpc about its centre giving a nominal 
resolution equivalent to 4096^ particles giving a mass resolution 
of rriDM ~ 2.95 x 10^/i^^Mq. Within this environment we iden- 
tified two main haloes, formally corresponding to the Milky Way 
and the Andromeda galaxy, whose main properties are listed in Ta- 
ble [T] together with their coiTesponding actual observational prop- 
erties. The virial mass of each halo is in units of IO^^A/q, while 
the virial radius and the distance between the two hosts, listed as 
D, are in Mpc. Both these quantities are based upon the definition 
M(< i?vir)/(47r/3J?^ij) = AvirPb whcrc p(, is the cosmological 
background density and Avir = 354 for the considered cosmol- 
ogy and redshift z = 0. The concentration is C2 = Rvir/r2, where 
r2 denotes the "scale radius" where the product p{r)r^ reaches its 
maximum value. The two axis ratios b/a and c/a are derived from 
the eigenvalues a > b > c of the moment of inertia tensor, and 
the vertical-to-planar axis ratio is reported for M3 1 . The a param- 
eter is the exponent corresponding to a scale-free host mass profile 
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Table 1. Main properties of the two lialoes (representing the MW and M3 1 
galaxy, respectively) considered in this work, and of their respective sub- 
haloes population. The virial mass of each halo is in units of IO^^A/q, 
while the virial radius and the distance between the two hosts, D, are in 
Mpc. We listed the observational in feiTed quantities o f M W and M3 1 , that 
refer t o the work of: (a) WIO, (b) iLaw et alj ll2009l). fc) iBaneriee & jog 
' 2008[) . (d ) iMcConnachie et alj i2005l), (e) iMated il998l). mllbataeTal 



20071) . f g) lMartin et alj fcOOSh . fh) lKarachentsev et al.l )2004h . 



property 


MW M31 


MW M31 




simulation 


observed 



Mvir 


1.674 


2.226 


2.7 iCS*"' 


1.5 ± 0.4('') 


-Rvir 


0.310 


0.340 


CSOO^"^) 


o.soot") 


C2 


11.7 


10.7 






b/a 
c/a 


0.937 
0.883 


0.978 
0.872 


0.83('') 
0.67('') 


0.4('=) 


a 


-0.034 


-0.052 






D 


0.782 


0.785 ± 0.025(''' 





1205 


1405 


24 


21 


Tout 


0.309 


0.340 


0.250 ± 0.003('=) 


0.270(-f) 


'^min 


0.018 


0.014 


0.023 ±0.002(9) 


0.005''^) 



AI{r) oc r^^", see Section IbTI for more details. The observation- 
ally derived masses are based on the work of WIO, and represent 
the estimates of each galaxy mass assuming a virial radius of 300 
kpc, using the observed anisotropy parameter /3 and including satel- 
lites' proper motions. 



2.2 The (Sub-)Halo Finding 

In order to identify haloes and subhaloes in our simulation we have 
run the MPI-l-OpenMP hybrid halo finder AHeO described in detail 
in|Knollmann & Knebe ( 2009). AHF is an improvement of the MHF 
halo finder jOill et al.l2004l) . which locates local overdensities in an 
adaptively smoothed density field as prospective halo centers. We 
would like to stress that our halo finding algorithm automatically 
identifies haloes, sub-haloes, sub-subhaloes, etc. and i t can reliably 
recov er substructures containing at least 30 particles (iKnebe et al.l 
l201lh . For more details on the mode of operation and actual func- 
tio nality we though refe r the re ader to the code description paper 
bv lKnollmann& Knebd t200Sl) , while an in- depth comparison to 
other halo finding techniques can be found in iKnebe et aljj201lh 
and Onions et al., in preparation. A complete summary of the char- 
acteristic of the subhaloes population of the two main haloes, MW 
and M31, is shown in Table [T] together with a comparison of the 
properties of their observed satellite galaxies. The rout and rin are 
the radius of the outermost and innermost tracer, respectively, in 
Mpc (in the case of M31 we listed the projected distances). The 
quantity A^sat represents the number of simulated subhaloes (or ob- 
served satellite galaxies) within 0.3 Mpc from each host center. 



(WIO) where the hierarchical structure formation model supports 
the notion that the density profile of dar k matter haloes follow s the 
functional form originally proposed by iNavarro et al.l (Il996h . i.e. 
the s o-called NFW pro file. 

IXue et al. I l l2008l) constrained the mass distribution of the 
MW's dark matter halo by analyzing the kinematic of thousands 
of blue horizontal-branch halo stars, finding a profile that is con- 
sistent with a combinat ion of a fixed disk and bulge model with a 
NFW dark matter halo. ISeigar et al.l fcOOSh have derived new mass 
models for M31, and found that while a NFW and an adiabatically 
contracted NFW profiles can both produce reasonable fits to the 
observed rotation curve of M3 1 , the pure NFW model requires a 
halo concentration too high with respect to the range predicted by 
the ACDM cosmology, and is therefore disfavoured. Thus, it is still 
debatable whether the Milky Way and Andromeda galaxy haloes 
actually follow a NFW profile. 

In this Section we briefly introduce the scale-free mass esti- 
mators, which are directly taken from WIO: we refer the reader to 
their work for a derivation of the respective formulae. 



3.1 Theory of mass estimators 

Here we present the four relevant formulae and the three parameters 
each formula depends on: we see that the mass estimator takes dif- 
ferent forms according to the available informations from the tracer 
populations. 



3.1.1 Full Information Estimator (FIE) 

In the optimum case that the full six-dimensional phase-space in- 
formation is accessible, the mass estimator can be written as: 



M(< 



GN ^ 



2 o 



with G = — TT^r „ut 



(1) 



(2) 



3-2/3 

Where v and r are the velocity and distance of each individual 
tracer particle, rout represents the radius of the outermost tracer, 
and G is the gravitational constant. The dimensional constant C 
is constructed out of three additional parameters determined by 
the host potential (a), the tracer's radial distribution (7), and the 
tracer's velocity anisotropy (/3), more details in Section l3. 1 .5 I where 
these parameters are algebraically defined. Note that we can only 
estimate the halo mass contained within the outer radius rout set by 
the distance to the farthest tracer. The mass is then constructed as 
an average of v^r°' over the total number of tracer objects, A^;i.acor- 
We will refer to Eq. ([T) as the Full Information Estimator or simply 
FIE. 



3 SCALE-FREE MASS ESTIMATORS 

Even though the mass estimators are derived under the assump- 
tion that the respective distributions are scale-free, they have nev- 
ertheless been successfully applied to the observed MW and M3 1 



3. 1.2 Radial Information Estimator (RIE) 

In the case that only the radial velocity, with respect to the center 
of the host galaxy, and the individual distances of the tracer popu- 
lation are known, Vr and r respectively, a different definition of the 
constant C must be used: 



AMIGA halo finder, to be downloaded freely from 

http : / /www . popia . f t . uam . es/AMIGA 



C 1 

M(< rout) = ""^'i 



(3) 
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3.1.5 The parameters a, (3, and 7 



with C = (a + 7 - 2/3)r^-t" 



(4) 



We shall call this the Radial Information Estimator, RIE: this case 
applies to our own Milky Way. Since we do not have the proper 
motion of all of it s satellites, but just of 9 of them (see for instance 
lMetzetalJ l [200^ ). it is safer to assume the RIE. It must be noticed 
that in absence of proper motion Vr may be calculated from vios by 
using the statistical correction: 



> 



(5) 



where ij> is the angle between the vector from the galactic centre to 
the satellite and the vector from the sun to the satellite. As we can 
see, this coiTection further depends from the anisotropy parameter 
/3. We will come back to the proper placement of the observer and 
the relevance of this correction, respectively, later on. 



3.1.3 Line-of-Sight Information Estimator (LIE) 

When using only projected line-of-sight velocities vios and actual 
distances r for the tracer population, the mass estimator referred 
to as the Line-of-sight Information Estimator, or LIE, and may be 
written as: 



The ever present constant C is composed of three parameters, de- 
scribing the host potential as well as particulars of the tracer pop- 
ulation, under the assumption that they both can be sufficiently de- 
scribed by scale-free models. We further assume spherical symme- 
try for our tracer population. 

The a parameter corresponds to a scale-free gravity field, 
which is equivalent to a host mass profile of the form: 



M{r)(xr°' (11) 
or, equivalently, to a mass density that scales as p oc r"'-""'"^-'. 

The /3 parameter is th e Binney velocity anisotropy parameter 
teinnev & Tremain3ll987l) . defined as: 



2cr2 



(12) 



in which a? and are the tangential and the radial velocity 
dispersions of the tracer objects. (3 provides information about the 
orbital distribution of our tracer population. 

Lastly, the 7 parameter represents the exponent of the power 
law describing the radial number density distribution n(r) of the 
tracer population: 



C 1 

M{< r-out) = 7^-7^ ^ vfos.iU, (6) 



. 3(a + 7 - 2/3) 

(7) 

This estimator must be used, for example, when calculating the 
mass of the Andromeda galaxy. 



3.1.4 Projected Information Estimator (PIE) 

In the worst case scenario in which the only data available are both 
projected distances R and line-of-sight velocities vios for the tracer 
population, the corresponding estimator is: 

M{<r^^,)^-- J2 ^fos,^R?, (8) 



with C = r^ut 

J-a,l3 

where 



n(r) oc r ^. (13) 

These three parameters are fundamental in describing the ge- 
ometry of the system and, together with the kinematical informa- 
tion of the tracers, allow us to compute an accurate estimator of the 
mass of a host halo. It is thus absolutely essential that they are de- 
termined with the highest possible accuracy. In reality however, this 
is not always possible: we are often forced to make assumptions re- 
garding the form of the underlying host potential. Moreover, the 
number of the known satellites of both Milky Way and Andromeda 
is only ~ 25, making the determination of the 7 parameter rela- 
tively inaccurate. In addition, in the Milky Way's case, only 7 of 
these objects have accurately measured proper motions: with such 
a small sample the velocity anisotropy (3 is widely unconstrained 
by data. In the next Section [T2l we will present the dependence of 
the mass estimator on each of these parameters, computing - for a 
specific case - the error introduced by uncertainty in a, (3, and 7, 
respectively. 



3.2 Dependence of the mass estimators on the parameters 

Given the inherent inaccuracy in determining the three model pa- 
rameters, we would like to gauge the sensitivity of the mass estima- 
tors to their uncertainties, considering an adequate set of subhaloes 
covering a radial range out to rout- 



ic.,f^= 4r(| + 5) + + (10) 

and r{x) is the gamma function. We will refer to this last equation 
as the Projected Information Estimator, PIE. 



3.2.1 J dependence 

We aim to study the dependence of the mass estimator on the pa- 
rameter 7, which represents the exponent describing the radial dis- 
tribution of the satellites population. We therefore calculated the 
relative variation of the estimated mass per variation in 7: 
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AM 



(14) 

MA-y a + 7 -2/3 ^ ' 

We note that Eq. l ll4b is valid for the four cases of FIE, RIE, LIE and 
PIE, being independent from the radial distribution of satellites; it 
provides a tool for calculating the expected uncertainty in the mass 
determination given the expected errors in 7. 

Assuming an isotropic distribution of orbits, i.e. /3 = 0, we 
focus on the real case scenario in which the 7 parameter is 7 = 2 
(as found in WIO for the observed satellites of MW and M31): 
allowing for an uncertainty of A7/7 ~ 25%, and recalling the 
typical value for a that is around 0.55 for a NEW host (e.g., WIO), 
we see that the error in the estimated mass is as high as AM/M ~ 
20%. This error will be even larger when considering smaller value 
of 7 and a, as well as for /3 > 0. 

3.2.2 P dependence 

Regarding the changes in the mass estimation due to the anisotropy 
parameter, we recall that this parameter is directly obtained from 
the velocities of the tracer population, computing the tangential and 
the radial velocity dispersion of each subhalo (as opposed to the a 
and 7 parameter, which are derived by assuming a power law dis- 
tribution). The P average value has been found to be /3 ~ —0.3 and 
P ~ —0.02 in our CLUES simulation, for the full set of subhaloes 
of the MW and M31, respectively. While these values of P indi- 
cate that we are close to the isotropic regime, i.e. P = Q, the MW 
anisotropy parameter slightly favors tangential orbits, in agreement 
with the measured proper motions of the known MW satellites. It 
is thus essential to understand how the variation in the /? parameter 
affects the determination of the host mass. For the RIE estimator 
the corresponding equation reads 



/ AM \ 



\MAp)^j^ Q + 7-2/3 
whereas for the FIE and LIE cases it is 



/ AM \ 
\MAp)^ 



F,E,L,E 3-2/3 a + 7-2/3 
finally, for the PIE scenario 



/AMX 
\MAp , 



a + 2 



(15) 



(16) 



(17) 



„g a + 3-/3(a + 2) Q + 7-2/J' 

In Fig. [T] we present the absolute value of the fractional mass 
variation as a function of the /3 value for the FIE and LIE cases 
(left panel) and for the RIE estimator (right panel). We do not plot 
the mass changes in the PIE case, as it is practically identical to 
the FIE and LIE ones. As in the previous section, considering the 
general case of having a NEW halo, with values of /3 close to zero 
and 7 = 2, which is the usual case for the hosts considered here 
and elsewhere (e.g., WIO), we find that the error due to variations 
of A/3 = ±1 for the FIE, LIE and PIE estimat ors is actually quite 
low and is below 10% for a — 0.55. Moreover. Ie vans et al. I( l201ll) 
found that for much of the radial regime covered by the tracer 
population, any variation of the anisotropy parameter within its 
physical range leads to the same estimator in the case of the PIE 
scenario. Thus, in the case of an external galaxy whose dark matter 
halo follows a NFW profile with a = 0.55 and 2 < 7 < 3, we 
can assume to have a minor error due to /3: the major uncertainty 



in the mass estimation comes from the assumption made on the 
a and 7 parameters. This last statement is valid for the FIE, LIE 
and the PIE estimator: it does not matter if we have real satellites 
distances or projected ones, the biggest error on the mass does not 
come from the anisotropy parameter. 

The situation is however, completely different for the MW 
galaxy, for which the RIE formula holds, i.e. we have radial in- 
formation on the satellite velocities. In this case, a variation of 
A/3 = ±1 could cause an error in the mass estimation of around 
80% if we have 7 = 2 and a = 0.55. Therefore, the /3 parame- 
ter is unfortunately the greatest concern in the calculus of the mass 
of our own Galaxy. Please note that if we knew the three dimen- 
sional velocities of the MW satellites as opposed to only the radial 
ones, we would be dealing with Eq. l ll6t . thus being in the regime 
in which the connect evaluation of the /3 parameter will only have a 
subordinate influence. 

We close by remarking that this discussion perfectly agrees 
with the previous study of the influence of the /3 parameter on the 
mass estimation as presented in WIO. 



3.2.3 a dependence 

Finally, we computed the amount or error introduced by uncertain- 
ties in the a parameter, which is directly connected to the potential 
of the host halo. The fractional variation of the estimated mass, for 
the FIE estimator, is; 



AM 
MAa 



1 



Q + 7 - 2/3 



ln(ro 



+ 



E. u^r" 
% % I 



(18) 

where the summation, as usual, is performed over the total number 
of tracers A^traccr and r-out is the radius of the outermost subhalo. 
Eq. 1 II8I 1 is formally identical for the RIE and LIE case as well, 
after substituting the full velocity v with the radial velocity Vr or 
the line-of-sight vios one, respectively. 

When dealing with the PIE scenario, instead, the error can be 
calculated through the following equation : 



AM 
MAa 



1 



a + 7 - 2/3 



In(rout) + h 



+ 



*(f + |) *(f + l) 



1-/3 



: + 3 - /3(q + 2) 



where '^(x) is the digamma function, defined as the derivative of 
the logarithm of the T(x) function. 

Unlike the other cases, we can not give a generalized estima- 
tion of the error introduced by the a parameter, it being dependent 
on the radial distribution of the satellites population: this uncer- 
tainty varies for every specific scenario and needs to be calculated 
individually. 



4 APPLICATION TO THE CLUES SIMULATION 

We now move to the application of the scale-free mass estimators to 
a situation as close as possible to our Local Group. To this extent we 
use the CLUES simulation introduced in Section |2] While we are 
certain that the scal e-free a pproximation le ads to credible result s 
as shown by WIO, HvanTet al. (2011)) and iDeason et al.l ( 1201 ih . 
it remains to be seen whether the uniqueness of the Local Group 
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Figure 1. Relative variation per unit parameter cliange of the mass estimated as a function of tlie /3 parameter in the case of FIE/LIE (left panel) and RIE (right 
panel) estimators. The fractional error in the mass estimation is larger for the RIE estimator, used for MW, than for the FIE/LIE estimators, applicable to M3 1 . 
See full text for more details. 



with its binary host sys tem and particular form ation history involv- 
ing preferential inf all jLibeskind et al.ll201lh . renegade satellites 
jKnebe et alj|201ll) and anisotropically distributed subhaloes (not 
explicitly shown here) will effect the mass estimate. Moreover, we 
would like to gauge the accuracy of these mass estimators when 
fewer tracers are used, as in the real LG. 

From now on we will refer to the case of an observer that is 
placed at the center of our Galaxy and looking towards the MW's 
satellites or to the nearby M31's ones. The choice to put the ob- 
server in the galactic center instead that at the solar radius may 
affect the determination of Vr from vias-, given the fact that the ra- 
dial velocity should be computed with respect to the sun. However, 
this is practically identical to the radial velocity with respect to the 
galactic center for distant tracers, for which sin ~ 0, which is 
the case for our subhaloes. Moreover, the anisotropy parameter /3, 
which appears in the correction factor of Eq. (|5), has always been 
found to be very close to zero in our simulations (as reported in 
Section [3.2.2| l. Nevertheless, when applying the different mass es- 
timators we also used the correction factor given by Eq. (|5), placing 
the observer on a sphere of radius 8 kpc from the galactic center, 
and we verified that the affect of this correction is at the < 0.5% 
level. We will thus refer, through the paper, to the case of an ob- 
server placed in the galactic center. 



assuming that the number density n{r-) follows Eq. ( 1131 ). 

The velocity anisotropy parameter /?, as defined in Eq. J12l (, 
is obtained by first calculating the radial velocity dispersion of the 
subhaloes, projecting their velocities along the radial axis, then by 
computing the tangential component of a through the relation 



/ 2 , 2 

= (o-j. 4- ay 



2\ 2 



(21) 



While P and 7 can be directly computed in the FIE and RIE cases, 
in order to calculate them in the LIE scenario we first need to 
derive the line-of-sight component of the velocity vectors of the 
subhaloes. The line-of-sight velocity depends on the viewing angle 
of the host which is unknown in our simulations. We thus randomly 
rotate each host and its subhaloes Nrot = 5 times, taking the 
mean of all these resulting line-of-sight velocity to compute /3. We 
perform a small number of rotations of the whole system since 
otherwise, by averaging over a higher number of rotations, we 
converge to the FIE case. The same methodology has been applied 
to the PIE case where we additionally had to project the distances 
of the tracers objects into the observers plane in order to obtain the 
7 parameter. 



4.1 Obtaining the parameters a, /3 and 7 

In order to apply the mass estimator method to our simulated galax- 
ies, we need to calculate the three unknowns a, /3 and 7 that appear 
in Eq. GJ,!©,© and Eq. {S}. 



4.1.1 The satellite parameters: 7 and P 

The 7 parameter is simply obtained by fitting the radial number 
distribution A'^(< r) of each host's subhaloes to the functional form 

iV(< r) oc r^""^, (20) 



4.1.2 The host halo parameter: a 

To get the value of a, we must recall that since our haloes are not 
scale-free but rather follow a NFW profile dNavarro et al.ll99^ , the 
applicability of a power-law is limited. While for a pure scale-free 
model it is irrelevant whether we fit the gravitational potential, the 
density or the mass profile of the host halo (see WIO), it will most 
certainly lead to differences when the scale-invariance is broken. 
Recall that for a scale-free model: 

(f>{r) oc <^ p{r) oc r^°^^ <^ M{r) oc r^^" . (22) 

For a NFW object however, we must identify which quantity is the 
most suitable to be fitted, and we decided to use the cumulative 
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Table 2. Value of the a pai'ameter and its fractional error Aa/a obtained 
by fitting the numerical mass profiles of the MW and M31 over different 
radial ranges. 





MW 




M31 




Radial Range 


a 


Aa 

a 


a 


Aa 

a 


[0, l]«vii- 


-0.034 


70% 


-0.052 


60% 


[0.4, l]R,i, 


0.302 


8% 


0.266 


9% 


[0.8, 


0.398 


3% 


0.402 


7% 



mass at the i?vir = 309 kpc within a 5% of uncertainty and the 
M31 mass at the i?vir = 340 kpc within a 3% of error, respectively 
(FIE estimator). However, we decided not to show these results and 
rather focus on more interesting and practical situations where the 
number of tracer objects is limited and agrees better with the actual 
observed Local Group. We must note however, that part of our ini- 
tial questions has been already answered by this exercise: the scale- 
free mass estimators are even applicable to a system of host haloes 
such as the (observed) Local Group for which they were originally 
designed. 



mass profile M (r) since this is the least noisy from a numerical 
point of view. 

Furthermore, since our halo does not follow a scale-free pro- 
file (either in mass or in potential), the actual value of a depends 
on the radial range used to fit it, i.e. a{r) 7^ const. 
We thus provide, in Table|2] the numerically fitted values of a, ob- 
tained by fitting the total mass halo profile in different radial ranges, 
specifically in the total range, in the outermost one, and in the inter- 
mediate range, together with their fractional relative errors, where 
we indicate with Aa the la error on a as found from the fitting 
routine. 

As in the previous case of the 7 parameter, we used a Poisso- 
nian weight (l/Naub) to associate errors to the data during the fit: 
as expected, the smallest relative error is obtained in the outermost 
radial range, confirming that in this regime the host density profile 
is best approximated as being of scale-free nature. 

We obtained for the MW and M31 in the total radial range 
a value of a = —0.034 and a — —0.052 respectively, as listed 
already in Table[T] while we can observe how the a value increases 
when we move to the outer part of the halo, as expected if the halo 
is following a NFW profile, since it gets steeper towards the outer 
part of the distribution. 

While using the numerical mass profile given by the simula- 
tion data is actually a self-consistent way to obtain q, we note that 
an observer would require a mass model to actually determine the 
Q parameter to be used with the mass estimators. Since an observer 
does not have any a-priori knowledge of the radial mass distribu- 
tion (or potential) of the host halo, an analytical profile must be as- 
sumed. Note that WIO showed that for an object following a NFW 
profile the typical value of a is ~ 0.55, based upon fitting a NFW 
potential in the range [10, 300] kpc to a power-law (p oc and 
assuming to have hosts with concentration between c = 18 and 
c = 8. Given the uncertainty on the actual density profile of the 
real hosts, we decided to allow for the estimates of a in two differ- 
ent ways: 

(i) using the values derived by fitting our numerical profile at 
different radial ranges, or 

(ii) using the relation a = 7 — 2, which holds true if the sub- 
haloes are tracking the total gravitating mass of the hosts. 

4.2 Results for the simulated MW and M31 

The application of the scale-free mass estimator to the (observa- 
tionally) unrealistic scenario in which we have A'^ ~ 1000 tracers, 
as found in our simulated haloes, gives excellent results for all the 
estimators FIE, RIE, LIE and PIE. Using the error formulae listed 
in Section [T2l and allowing a maximum error on the calculation of 
the parameters a, /3 and 7 of about ~ 20%, we obtained the MW 



4.2.1 Matching the number of the obserx'ed satellites 

As shown in Table [T] the total number of subhaloes found within 
300 kpc in our simulations substantially differs from the number 
of observed satellites galaxies of the Milky Way and Andromeda 
within the same ra dius (the well known miss ing satellites p roblem , 
first addressed in iKlvpin et all h999l) and iMoore et all ( Il999l) ). 
Thus, we would like to calculate the accuracy of the mass estima- 
tors when the number of tracers is comparable to the real one, i.e. 
A'^ ~ 30 (we explicitly chose this number to be able to have a di- 
rect comparison with the WIO results, see for example their Fig.l). 
Further, the real case scenario is the one for which we have the ra- 
dial velocities of the MW satellites and the line-of-sight velocities 
of the M31 tracers: in the forthcoming analysis we will thus only 
use the RIE estimator for the Milky Way and the LIE one for the 
Andromeda galaxy. 

From the total set of subhaloes we randomly selected A" = 30 
objects that covered the total range within rout < 300 kpc and 
computed their velocity anisotropy and their radial distribution, 
thus obtaining the /? and 7 coefficients. For this particular exercise, 
the a parameter was numerically evaluated using the three different 
radial ranges of the host mass profile listed in Table|2](ignoring the 
option to evaluate it as a = 7 — 2 for the moment). For each of these 
values of a we performed 1000 random realization, we applied the 
scale-free mass estimator and we calculated the distribution of the 
ratio of the estimated over the actual mass, i.e. A/est/A/truc; those 
distributions have then been fitted by a Gaussian curve eventually 
leading to the best-fit parameter /i and its standard deviation a. 

The results of these tests (for the a value evaluated from the 
total radial range, i.e. first line of Table|2) are summarized in Fig.|2] 
for the MW and Fig. [3] for M31 where we plot in the left panels 
the distributions of Mcst /Mtmc for the FIE mass estimators and in 
the right panels the RIE (MW) and LIE (M31), respectively. The 
legends of each panel further list the three parameters a, 7, and /? 
relevant for the respective mass estimator (where /3 and 7 repre- 
sent the average value over the total 1000 realizations) alongside 
the peak and standard deviation of the best fit Gaussian. Note that 
the standard deviation is compatible with 1 /\/N7ub where Ng^fy is 
the number of used tracers, and it increases when only radial veloc- 
ities (or line-of-sight ones) are used. Remarkably, the mean of the 
distribution stays always very close to fj, = 1.0: the mass estima- 
tors are thus unbiased with respect to the number of used objects. 
We repeated the above mentioned analysis for the other values of 
a listed in Table|2] and found practically indistinguishable results: 
the best-fit ^ and a values are given in in Table |3] 

In summary, we found that for both host systems the mass is 
always recovered within a few percent of error when restricting the 
analysis to 30 randomly selected subhaloes each. 
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Figure 2. Distribution of 1000 realizations of the estimated mass over the real one for 30 subhaloes of the simulated Milky Way. The average value of the 
parameters /? and 7 obtained in each realization is shown. The best-fit Gaussian is also plotted, and its mean ^ and standard deviation a are indicated. The left 
panel corresponds to the FIE estimator, the right panel to the RIE one. 
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Figure 3. Distribution of 1000 realizations of the estimated mass over the real one for 30 subhaloes of the simulated Andromeda galaxy. The average value of 
the parameters /3 and 7 obtained in each realization is shown. The best-fit Gaussian is also plotted, and its mean ^ and standard deviation a are indicated. The 
left panel corresponds to the FIE estimator, the right panel to the LIE one. 



Table 3. Mean value and standard deviation of the Gaussian distribution of 
the estimated mass over the true mass for the MW and M3 1 hosts, using 
the RIE and the LIE, respectively, for the three considered values of a (cf. 
Table|2j. The number of subhaloes has been limited to 30 random ones, and 
they have been used to compute /3 and 7 for each realization. The average 
values of these parameters over the total N = 1000 realizations are 7 = 
1.63±0.12and/3 = -0.307±0.061 for the MW, and 7 = 2.013±0.013 
and/3 = -0.006 ± 0.001 for M31. 





MW 






J\/31 




a 




a 


a 






-0.034 
0.302 
0.398 


1.057 
0.992 
0.975 


0.204 
0.177 
0.169 


-0.052 
0.266 
0.402 


1.060 
0.958 
0.931 


0.257 
0.167 
0.159 



4.2.2 Matching the radial number distribution of the obser\'ed 
.satellites 

While using 30 randomly chosen subhaloes leads to exceedingly 
well recovered host masses, we acknowledge that our model sub- 
haloes (for the MW) do not follow the same radial distribution as 
the observed ones (why this is the case is substance for yet an- 
other paper and shall not be addressed here): we list in Table |4] 
the distances to all present ly known MW satellites (taken from 
IWadepuhl & Springelll20r 11) alongside their masses and use this 
data to obtain the observed 7 ± A7 = 2.80 ± 0.08 by fitting the 
radial distribution to a power-law in Fig. |4] Please note that we 
only focus on the MW's subhaloes here, as in the case of M31 the 
7 = 2.013 coefficient is very similar to the one obtained from the 
observed satellites distribution (see WIO). 

From the total set of subhaloes in our numerical MW, we con- 
structed a subset of 30 tracers by selecting those objects that follow 
the radial distribution 7V(< r) oc r''"^'*. Further care was taken to 
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Figure 4. Radial distribution of the observed MW satellites and coiTespond- 
ing best fit in the range r < 300 kpc. 



verify that the randomly selected subhaloes always cover the (ob- 
servational) radial range up to ~ SOOkpc. While the 7 = 2.8 is 
fixed by construction the /3 has always been derived from this sub- 
set; for the a we first used, again, the three values listed in Ta- 
ble |2l and we found a notable bias in the Gaussian distribution of 
Mest/Aftruc, as high as the 80%: this choice of a does not pro- 
vide the expected host mass. Thus, we secondly decided to verify 
if the assumption that the tracers are tracking the total gravitating 
mass of the host can provide a better constrain on the value of a, 
i.e. using the relation a = 7 — 2. In this case, without making any 
fits to the numerical shape of the host profile, we actually found re- 
sults in excellent agreement with the true mass, as shown in the left 
panel of Fig. |5] In the right panel of the same figure we show the 
distribution obtained when yet another additional constraint was 
added, i.e. we selected only those subhaloes whose mass lies be- 
tween 5 ■ 10'' < M/Mq < 1 • 10®, in order to resemble the aver- 
age mass of the observed MW satellites (see Table |4j. Also in this 
case we can observe that the Gaussian is peaked very close to 1 .0, 
at/x = 1.016. 

We finally decided to also test and use the suggested value of 
a = 0.55 (WIO), but we actually obtained a Gaussian mean value 
for fi that is biased by approximately 30% towards large estimated 
masses. 

In summary, even when restricting the subhaloes to follow the 
same power-law as the observed satellites within the same mass 
range, the scale-free mass estimators are capable of recovering the 
true mass of our constrained MW and M31 if one chooses to use 
a — -y — 2 (being close to the isotropic regime, i.e. /? = and as 
far out as /3 = —0.5). 

4.2.3 Do we require a host mass profile or simply a = ^ — 2? 

The analysis in the previous subsection has shown that simply using 
a = 7 — 2 actually leads to excellent results for the scale-free mass 
estimators when applied to our constrained Local Group and a sub- 
halo population restricted to follow the observed one as closely as 
possible. But can this finding be generalized, at least with respects 
to our simulation? 

To verify that the assumption a = 7 — 2 holds true in general, 
we select the subhaloes of both MW and M3 1 in order to follow 
different radial distributions, according to A''(< r) oc r"^^^, where 
we allowed the 7 to vary between 1.5 and 3.0. The a coefficient 
was then calculated consequently, while the /3, as usual, came from 
the selected satellites velocity dispersions. In this way we selected 



Table 4. List of the MW satellites used in this work, con'esponding to those 
lying within 300kpc from the galactic center and with measured line-of- 
sight vel ocities. The Galacto cent ric distances D a re in kpc. The values ar e 
from: (a)lMartin et al.ll l2O08l) . ( b)lMatedlll998h (c)lBelokurov et allfcOOSh 
(d) Iva nden Bergh ( 199i). ( e^ Ifielokurov et ajj j2009l) . (f) ISimon & GelS 
{im'h . (g) lBekkil l2008i) . (h) lvan den Berghl j2OO0l) . 



Name 


D[kpc] 


Mass[10^MQ] 




66 ±3 


- 


Booll'-'^') 


42 ±8 


- 




101 ± 5 


13 




44 ±4 


1.2±0.4(^) 


CVnI<-°-'> 


218 ± 10 


27±4'-f'> 


CVnll'-'^'l 




2.4± l.lt-'') 


Draco^"^ 


76 ±5 


22 




138 ± 8 


68 




132 ± 12 


7.1 ± 2.6(^) 


LeoI<-'''> 


250 ± 30 


22 


Leoll'-''^ 


205 ± 12 


9.7 


LeoIV'-'^'l 


16011^ 


1.4± 1.5(^') 




180 




LMC'-'^') 


49 


10.000(9) 




24 ±2 


150''') 


SculptoA''') 


79 ±4 


6.4 


Segl'-'^'l 


23 ±2 




Segll'-'^') 


35 




Sextans^'''' 


86 ±4 


19 




58 


400(9) 


UMaI<-'^'l 


96.8 ±4 


15±4(^) 


UMaIl(°-) 


30 ±5 


4.9±2.2(-'') 




66 ±3 


23 




38 ±7 





A'^ = 30 subhaloes again for 1000 times and we obtained Gaussian 
distributions of the Mcst/A/truo quantity. We show the best-fit fi, 
with the standard deviation a as error bars, as a function of 7 in 
Fig.|6]for the MW (left panel, RIE only) and M31 (right panel, LIE 
only). The anisotropy parameter was always found to be close to 
/3 ~ 0, with a maximum variation between —0.3 < (3 < 0.1, in- 
dicating that we are in the isotropic regime. We would like to high- 
light that despite other choices of a may in principle be possible, 
as demonstrated in Section |4.2.1| the simple assumption a = 7 — 2 
provides always the best estimation for the host mass, whose as- 
sociated Gaussian distribution has mean value p closer to I and 
smaller standard deviation a. 

We see that the simple assumption, that the satellites track the 
total mass of the host halo, is sufficient to give an excellent es- 
timation of the parameter a to be used. The suggested value of 
a — 0.55, indicative of a NFW halo potential, is thus recommend- 
able when the satellite distribution follows a power law with expo- 
nent 7 ^ 2.5 : these valu es have been alread y successfully used in 
lOeason et alj j20Ill) and Ie vans et alj 1 I2OI ll) . This is of fundamen- 
tal importance for observations: without having any a-priori knowl- 
edge about the host halo density profile, we can simply use the 
value a = 7 — 2 once we have calculated 7 from the satellite dis- 
tribution. This condition has been verified in our constrained sim- 
ulations only, in which the anisotropy parameter is always /3 ~ 0: 
care should be taken when dealing with satellite populations whose 
P is highly anisotropic. 

In summary, we have shown that our two model hosts within 
the simulated constrained Local Group allow the application of 
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Figure 5. Distribution of the estimated RIE mass over the real one for 30 subhaloes of the Milky Way and 1000 realizations. The subhaloes have been selected 
following the power law with 7 = 2.8 (left panel) as well as additionally also constraining them to lie within the observed mass range (right panel). Note that 
in both cases a has been determined as a = 7 — 2. 



scale-free mass estimators to them. And, for as long as we are in 
the isotropic regime in which /? = 0, the simplifying assumption 
of a = 7 — 2 can be used. This alleviates the need to derive this 
parameter from a model of the host potential or mass profile. 

4.3 Exploring the influence of the MW and M31 on each 
other 

4.3.1 Subhaloes (not) facing the opposite host 

The fact that the MW and M3 1 hosts are close to each other, form- 
ing a binary galactic system, poses the question if the mass estima- 
tors work even in the situation in which we only consider satellites 
in between the two haloes. To shed light on this issue we begin by 
separating the MW halo into two hemispheres, defined as "facing" 
and "non-facing" M31 (we perform the same test for the M31 halo, 
too). Subhaloes are then grouped by the hemisphere they sit in. We 
remark that this is a purely spatial cut, to investigate if the proxim- 
ity of the companion host causes some bias on the mass estimator. 

The facing/non-facing subhaloes of the MW are selected ac- 
cording to a radial number distribution A'^(< r) oc r'^"^ with 
7 = 1.63, and consequently a = 7 — 2 = —0.37, while the 
M31 subhaloes are selected in order to follow the distribution with 
7 = 2.013 and a — 0.013. We chose these value to match the 
parameters found in Section 14.2.1 1 but verified that our results are 
not affected by the choice of the specific power law, as already 
expected (cf. Fig. [6}. Finally, we randomly selected A'^ = 30 sub- 
haloes from each of the two facing/non-facing subsets, repeating 
the analysis 1000 times and computing each time the parameter (3 
and the estimated mass. The anisotropy parameter, for both MW 
and M31 and in every subset of objects considered, has been found 
to be very close to again, lying between —0.35 < /3 < 0.05. We 
are hence in a situation to explore the influence of the two hosts 
on each other: with the /3 close to and the subhaloes selected to 
follow a fixed power law (thus without errors associated to the 3 
main parameters) we can affirm that any deviation in the mass esti- 
mation should now be attributed to the subhaloes facing or not the 
other host. 

In Table|5]we show the results of the mean value and standard 
deviation of the Gaussian distribution of the estimated mass over 



the true one for the MW (by applying the RIE) and for M3 1 (using 
the LIE). We show the /i obtained by using the facing subhaloes, 
the non-facing subhaloes and also the total set. We did this exercise 
for two different values for Tout, thus computing the host mass up 
to this outer radius: in this way we should be able to observe if the 
proximity or, on the contrary, the distance of the subhaloes popula- 
tion to the opposite host has some influence as well. For the MW 
host we can observe that the estimator performs equally well when 
using the facing or non-facing objects, for each of the two rout val- 
ues used. In the case of M31, instead, the non-facing population 
of subhaloes seems to give better results in the estimation of the 
host mass, while the facing objects lead to a Gaussian distribution 
whose mean value is slightly biased (p — 1.11) when we consid- 
ered the Tout = 221 kpc. Given the high uncertainties associated to 
this biased result however, we can conclude that each of the main 
galaxies does not influence the subhalo population of the other one. 



4.3.2 Renegade subhaloes 

As already discussed, we call renegade subhaloes those objects, 
in our simulations, that change their affiliation fro m one of the 
two p rominent hosts in the Local Group to the other dKnebe et all 
I2OIII) . We were able to identify 129 renegade objects, 57 of which 
belonging to M31 at z=0, and the remaining 72 belonging to the 
MW. 

We thus examine the effect that this population of renegade 
subhaloes may have on the mass estimators: while in the previous 
studies we used the full set of subhaloes, automatically including 
also the renegade ones, we now want to restrict the analysis to the 
renegade subhaloes only in order to estimate the mass of MW and 
M31. 

For each host we computed the anisotropy parameter and ra- 
dial distribution of their respective renegade satellites, and found 
Pmsi = 0.86 - Pmw = 0.72 and 7M31 = 2.15 ± 0.21 - 
7A/W = 2.01 ± 0.20. The choice of the host parameter a is made 
considering its value in the total [0, 1]-Rvir range, or using the re- 
lation a = 7 — 2 that we provided in the previous sections. The 
resulting estimated masses are shown in Table |6] in which we have 
used the FIE estimator in order to compare the effects of these rene- 
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Figure 6. Mean value of the best fit Gaussian distribution for different values of 7 and correspondingly = 7 — 2. The anisotropy pai'ameter lies between 
—0.3 < /9 < 0.1. The distribution is based upon 1000 realization. MW, using the RIE, left panel, M31, using the LIE, right panel. 



Table 5. Milky Way and Andromeda galaxy mass estimation (RIE for MW, LIE for M3 1 ) using the subset of subhaloes facing and non-facing the companion 
host. The mean value ^ and standard deviation a of the Gaussian distribution of the estimated mass over the true one are shown, obtained by selecting N=30 
subhaloes from the facing or not facing subsets, repeating 1000 realizations and evaluating the mass at two different rout. 



Host Tout 


All Subs 


Facing Subs 


Not Facing Subs 


'"out,MW 


/X ± (T 


/I ± (T 


/X ± CT 


309 kpc 
197 kpc 


0.97 ±0.12 
0.98 ±0.12 


1.05 ± 0.16 
1.01 ± 0.12 


0.85 ±0.24 
1.03 ±0.26 


''out.MSl 


/.i ± cr 


± (T 


/X ± CT 


340 kpc 
221 kpc 


1.02 ±0.19 
1.05 ±0.19 


1.08 ± 0.19 
1.11 ± 0.22 


0.92 ±0.15 
0.98 ±0.16 



gade subhaloes in the same way for both hosts. This time, because 
of the small number of objects considered, we do not perform mul- 
tiple realizations of the mass estimation, but only one: the errors 
associated with the mass are thus computed through the error prop- 
agation formula based on Eq. ( I14t - Eq. M9\ . in which we further 
assume that /3 is fixed, A7/7 ~ 10% as obtained by the fitting 
routine and Aa/a is listed in Table |2] We see that in the case of 
using the value of a from the total radial range, the mass estimator 
is biased for both hosts, with a large associated error. When us- 
ing the relation a = 7 — 2, instead, the mass of both Milky Way 
and Andromeda in recovered within a much smaller uncertainty. Is 
interesting to note how the relation between the host parameter a 
and the subhaloes distribution parameter 7 seems to hold true also 
in this case, in which the anisotropy parameter /3 is substantially 
far from being isotropic. However the lack of statistic in this case, 
having at our disposal just one realization of a small number of 
renegade objects, prevents us from generalizing the finding of Sec- 
tion |4j23] to this highly anisotropic case. The fact that /3 ~ 1 for 
the renegade subhaloes means that these objects are mainly mov- 
ing on radial orbits with respect to their hosts. We conclude that the 
computation of the host mass based upon a family of renegade sub- 
haloes gives results in agreement with the expected ones and hence 
these mass estimators will not be biased in case that renegade ob- 
jects also exist in reality. 



Table 6. Estimation of the MW and M3 1 mass based upon renegade sub- 
haloes only, FIE case. The parameter used are /3 = 0.72 and 7 = 2.01 for 
the MW, f} = 0.86 and 7 = 2.15 for M31. We remind that the value of a 
from the total radial range is a = —0.034 for the MW and a = —0.052 
forM31,asinTable|2] 



MW M31 


a 


{Mest ± AMest)/Mtrue 


[0,l]flvir 


0.92 ±0.33 0.75 ±0.40 


7-2 


1.00 ±0.21 0.97 ±0.16 



4.3.3 Unbound subhaloes 

For all previous results we did not test whether or not a subhalo 
is gravitationally bound to its host; we simply used a spatial cri- 
terion to determine its affiliation as this is how satellites are often 
defined in the observations. Now instead, we impose an additional 
constraint on our subhalo population: its velocity has to be lower 
than the local escape velocity Veac of the halo at that distance. Fol- 
lowing this criterion, we find that about the 3% of the subhaloes 
inside the virial radius of each host are unbound. As expected, most 
of them lie near by the virial radius, where the Vesc is lower and the 
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effects of the proximity of the other host are more important. We 
thus quantify the effects that unbound subhaloes have on the mass 
estimators. This is an interesting test, as it corresponds to the real 
case scenario in which the affiliation of a tracer object is not clear 
and could be erroneously included into the calculation of the mass 
of a host. 

We repeat our previous methodology by evaluating the MW 
and M31 mass 1000 times with a subset of A'^ = 30 subhaloes, 
this time including one, two or three unbound subhaloes. In order 
to ensure an unbound subhalo is included, we explicitly substitute 
in each realization, one, two or three of the 30 subhaloes with an 
unbound one. We then computed the /3 and 7 coefficients for each 
realization, and used the formula q = 7 — 2. We verify that the in- 
clusion of a single unbound subhalo leads to mass estimators which 
are slightly biased towards larger masses: we obtained, for both 
M31 and MW, a Gaussian peaked at ^ = 1.04 with a = 0.14. 
When including two unbound subhaloes, we found an higher devi- 
ation, with ^ = 1.12 and o = 0.13. Finally, forcing three unbound 
subhaloes to be included in each 30 subhalo subsample, we obtain 
a Gaussian peaked at /x = 1.17 and cr = 0.14. We performed the 
calculation using the FIE estimator but verify that our results are 
the same in the RIE and LIE case, giving similar values for the 
mean of the distribution /i, and increasing standard deviations a 
with respect to the FIE case. 

We note that the results presented in the previous sections are 
not significantly affected by the presence of unbound subhaloes, as 
in that case the probability that in a given realization of A'^ = 30 
randomly picked subhaloes one was unbound is due 
to the fact that unbound subshaloes make up just 3% of the full 
subhalo population. In this last test, instead, the probability that 
one object is unbound, over the A'^ = 30 subhaloes used in each 
realization, is P = 100%, because we deliberately replaced one 
random subhalo with an unbound one. Thus, we expect that the er- 
ror on the mean value of the Gaussian distributions in the previous 
analysis, caused by the possible inclusion of one unbound object, 
is 100/37 ~ 2.7 times lower than the error made here, when one 
subhalo is forced to be unbound. Looking at the fi = 1.04 obtained 
in this section, for the FIE case when we used a single unbound 
subhalo, we see that the 4% of deviation from the expected value 
will be reduced of a factor 2.7, giving negligible eiTors. We are fur- 
ther reassured by having performed the analysis of Section l4!2l also 
by explicitly neglecting the unbound objects, and we have observed 
no significant differences in the results already presented. 

To conclude, in this section we demonstrate that, being sure 
of having included unbound subhaloes, this inclusion causes an 
overestimate of the host mass, in agreement with the results of 
iDeason et all l l201lh . The more unbound objects we include into 
the mass estimator, the more biased the final mass is. Care should 
thus be taken when considering objects at the "edge" of a galaxy 
halo, as they may be not bound to it. 



* This probability can be computed using the hypergeometric distribution, 
which, in our case and for the M31 and MW, respectively, describes the 
probability to get one unbound subhalo within fc = 30 randomly drawn 
objects from a total subhaloes population of size A'^ = 1405(1205) in 
which the unbound objects are n = 45(36), thus the ~ 3% of the total. 



OC::)' 



■.kJ 



(23) 



5 SUMMARY & CONCLUSIONS 

We verified t he accuracy of the sc ale-free mass estimators recently 
proposed bv IWatkins et al.l ( |2010| . WIO) when applied to the two 
dominant Local Group host haloes, the MW and M31, by using 
Constrained Local UniversE Simulations (CLUES). These scale- 
free mass estimators assume that all the relevant informations about 
the enclosed mass of a halo are contained in the properties of its 
satellites (or any other tracer population), such as distances and ve- 
locities, which can both be given as either projected or full 3D 
data. The importance of such estimators resides in the fact that the 
full six-dimensional phase-space information of all celestial bod- 
ies down to the very faint magnitude G ~ 20 mag will soon be 
available thanks to the upcoming GAIA mission. What makes the 
usage of these mass estimators so appealing? After three years of 
operation the Sloan Digital Sky Survey II (SDSS-I0) discovered 
eight new dwarf galaxies, seven of them orbiting our Galaxy. The 
SDSS, which covered more than a quarter of the sky, essentially 
doubled the known number of MW satellite galaxies, helping close 
the gap between the observed number of dwarf satellites and the- 
oretical predictions. During its projected five-year mission, GAIA 
will scan the entire 41253 square degrees of the sky, obtaining as- 
trometric parameters (angular position, pr oper motion, an d paral- 
lax) for roughly one billion stars. Recentlv. lAn et al.](l2012h further 
investigated the benefits that the use of all this new proper mo- 
tion data will introduce in the application of mass estimators. It is 
thus imperative to develop and test against simulations the mass 
estimator based entirely upon tracer objects, such as s atellite galax- 
ies. Th is is sue has already been partially addressed in lDeason et al.l 
( l201lh and lEvans eTal] ( 1201 ih . using the GIMIC suite of simula- 
tions, from which they selected a set of galaxies that resemble the 
Milky Way. 

In this work we tested the scale-free mass estimators against 
Constrained Simulation of the Local Group, in which observational 
data of the nearby Universe is used as constraints on the initial 
conditions. These constrained simulations provide us with a unique 
Local Group seated in the correct environment, as opposed to un- 
constrained cosmological simulations, to make a direct comparison 
between numerical results and observations: verifying the goodness 
of the WIO mass estimators in our simulated LG should therefore 
been seen as complementary to the already published work on their 
credibility and as an extension to a system resembling as closely as 
possible the real Local Group. 

Our motivation is driven by the fact that the Local Group 
likely is a rather speci al (binary) system of galaxies featur- 
ing backsplash galaxies l lKnebe et alj 1201 il), renegade satellites 
(Knebe et al. 2 01 1| ) and preferential infall: Libeskind et alj 1 I2OI ih 
have studied the simulated MW and M31 galaxies in the CLUES 
framework and have found a clear evidence for the anisotropic in- 
fall of subhaloes onto the ir respective hosts. This result has been 
recently conoborated by iKeller et"aL 1 ( |2012|) . who examined the 
spatial distribution of the MW young halo globular clusters find- 
ing that they are anisotropically spatially distributed, sharing the 
same accreted origin as that of the MW's satellite galaxies. Our 
simulations also show the typical anisotropy in the distribution of 
subhaloes population, compat i ble with the obse r ved cl assical MW 
satellites jkroupa et allbOOSl : iMetz et alj|200l |2008|). and as al- 
ready found in other cosmo l ogical simulations ( cf. iKnebe et al.l 
|2004 iLibeskind et"ZI l2005l : IZentner et all l2005h . Furthermore, 
when comparing constrained against un-constrained simulations - 
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only 1-3% of Local Group (candidate s) share similar formation 
properties jporero-Romero et al.lbou]) . Thus, it is clear that our 
Local Group is a very special object in the Universe. 

We first studied the sensitivity of the mass estimators with re- 
spect to their main parameters: a, which describes the host halo 
scale-free gravity field, /3, which corresponds to the satellites' ve- 
locity anisotropy, and 7, representing the exponent of the power law 
describing the number density of the tracer population. We found 
that for an external galaxy, such as M31, for which only line-of- 
sight informations are available, the greatest error comes from the 
uncertainty of a and 7, whereas the mass variations stemming from 
the anisotropy parameter /3 are almost negligible in the interesting 
physical range. On the other hand, the greatest concern in the esti- 
mation of the mass of our Galaxy comes from the /3 parameter, as 
pointed out by WIO: without precise information about the satel- 
lites' proper motions, the error introduced by using their radial ve- 
locities is sensibly high. Hopefully, future surveys (e.g. the GAIA 
mission) would be able to measure such proper motions, signifi- 
cantly improving the quality of these mass estimators. 

We then applied the relevant mass estimators to the MW and 
M31 Local Group system as found in our constrained simulation. 
We found that all the estimators (FIE and RIE for the MW, FIE 
and LIE for M31) provide an unbiased results, with the mean of 
the Mest/Mtrue distribution close to /i = 1.0 and its standard 
deviation scaling with Nsub, even when a small {N ~ 30) 
number of tracer objects are used, resembling the real case scenario 
of the known satellites galaxies. When selecting the subhaloes in 
order to follow a specific radial number distribution 7V(< r) oc 

' with different 7, we found that, in the limit of the isotropic 
regime (i.e. /3 ~ and as far as /? = —0.5 in our simulations) 
the assumption that the subhaloes are tracking the total mass of the 
host (i.e. a = 7 — 2) is sufficient to get a very good estimate of the 
host mass. 

We thus also investigated how the mass estimators work when 
using subhaloes that are closer or further away from the neighbor- 
ing host, by restricting the analysis to the facing and non-facing 
hemispheres and calculating the mass at different values for rout: 
we observe that the two dominant hosts of the Local Group do not 
appear to influence its subhalo populations - at least not when the 
applicability of the scale-free mass estimators is concerned. 

Finally, we explored the possibility that using renegade sub- 
haloes, i.e. subhaloes that change their affiliation from one of the 
two hosts to the other, in the estimation of the mass may cause a 
bias: we do not find evidence for this, on the contrary, we observe 
that the mass of both MW and M31 is recovered with a few per- 
cent of error when the assumption a = 7 — 2 is made, even if the 
anisotropy parameter in this case is /3 ~ 0.7, indicating that these 
objects are moving on radial orbits. 

On the other hand, the inclusion of unbound objects, mainly 
found near the virial radius of each host, is able to cause an over- 
estimate of the host mass, as high as the 20% when 3 unbound 
subhaloes are used out of a total of 30 objects. In this regard, care 
should be taken when dealing with tracer populations whose af- 
filiation to each of the two host is not clear. As long as bound- 
ness is verified however, the unique subhaloes population in our 
simulations, sh owing anisotropy in th eir spatial distribution, pref- 
erential infall jlTbeskind et al.ll201 lb and even renegade objects 
jKnebe et alj|201ll ), does not prevent us from always recovering 
a good estimation of the host mass. 

Hence, the most important finding of this work is that satellite 
galaxies are well suited to "weigh" the MW's halo. Even with a 
small sample of just two dozen objects and despite anisotropic ac- 



cretion, an anisotropic spatial distribution, different masses, sizes 
and histories, subhalo kinematics are dominated by the host poten- 
tial, making satellite galaxies well suited for the problem at hand. 
We therefore conclude that the application of the scale-free mass 
estimators to either of the two dominant hosts of the Local Group 
provide credible results, it therefore appears safe to use it for the 
Local Group as already done by WIO. To get a good estimation of 
a host mass, in the case of having the anisotropy parameter /3 ~ 0, 
we recommend the use of the parameter a = 7 — 2 where the 7 
directly come from the observation of a satellite population. 

In the future, sensitive surveys and space based telescopic mis- 
sions will most likely improve both the census of satellite galaxies 
as well as our understanding of their proper motions: these new 
data will enhance the mass estimators making them more accurate 
than they are today. 
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